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The Tethered Particle Motion (TPM) technique informs about conformational changes of DNA 
molecules, e.g. upon looping or interaction with proteins, by tracking the Brownian motion of a 
particle probe tethered to a surface by a single DNA molecule and detecting changes of its amplitude 
of movement. We discuss in this context the time resolution of TPM, which strongly depends 
on the particle-DNA complex relaxation time, i.e. the characteristic time it takes to explore its 
configuration space by diffusion. By comparing theory, simulations and experiments, we propose a 
calibration of TPM at the dynamical level: we analyze how the relaxation time grows with both 
DNA contour length (from 401 to 2080 base pairs) and particle radius (from 20 to 150 nm). Notably 
we demonstrate that, for a particle of radius 20 nm or less, the hydrodynamic friction induced by the 
particle and the surface does not significantly slow down the DNA. This enables us to determine the 
optimal time resolution of TPM in distinct experimental contexts which can be as short as 20 ms. 



INTRODUCTION 

Biophysical techniques at the single molecule level have 
become an integral part of the available toolbox to in- 
vestigate biomolecular machineries. The rapid develop- 
ment of experimental techniques for the exploration of 
conformations and dynamics of single DNA molecules 
emphasizes the need for suitable theoretical tools to in- 
terpret the large amount of data collected. Among the 
many techniques used in biology or biophysics labora- 
tories. Tethered Particle Motion (TPM) [IUl3] is very 
promising because it explores the equilibrium statisti- 
cal mechanics of the biopolymer in the absence of ex- 
ternal force, by contrast to optical or magnetic tweez- 
ers experiments [14]. One end of a DNA molecule is 
immobilized on a glass surface, and the other end is 
attached to a particle, the diameter of which ranges 
from a few tens to several hundreds of nanometers (fig- 
ure [I]) . By measuring DNA end-to-end distance (or "ef- 
fective length"), tracking by video-microscopy the par- 
ticle trajectory informs about DNA conformations in 
real time. Hence TPM also gives access to dynami- 
cal properties. Getting insights into the dynamics of 
biomolecular events is of great significance. Should its 
time resolution be sufficient, TPM has the capability 
to get access to the dynamics of DNA conformational 
changes, such as looping/unlooping [4l[6l I13j. curvature 
variations due to protein binding/unbinding [HHS], hy- 
bridization/dehybridization [TB], or changes induced by 
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enzyme processing [H, ^ . Such conformational changes 
are detected through variations of the particle amplitude 
of movement, which is calculated on sliding time intervals 
of a given duration, Tav. On the one hand, Tav must be 
long enough to have a good estimate of the amplitude of 
movement, and thus to discriminate different amplitudes 
associated with different DNA conformations. On the 
other hand, Tav sets the TPM time resolution because 
events shorter than Tav a-i'e smeared out and thus cannot 
be detected. Therefore Tav must be optimally chosen. 
We shall see below that its optimal value is proportional 
to the relaxation time of the particle-DNA complex, i.e. 
the characteristic time the complex takes to explore its 
configuration space by diffusion. Thus the knowledge 
of the relaxation time is a prerequisite to estimate the 
TPM time resolution. The experimental conditions, es- 
pecially the proximity of the surface and the attachment 
of the particle, are likely to perturb the polymer and to 
increase the relaxation time. Therefore knowing the re- 
laxation time in function of both the DNA contour length 
and the particle radius is of primary importance to an- 
ticipate TPM time resolution capabilities. 

To our knowledge, dynamical consequences of the 
setup geometry, in particular of the attached particle, 
have not been quantified extensively yet. Calibration of 
TPM experiments, at the dynamical level, remains to 
be performed. We focus on the DNA relaxation time, 
technically defined as the characteristic time of the slow- 
est elastic mode only. We analyze how it scales with 
particle radius and DNA length, by using small parti- 
cles with radii down to 20 nm, for DNA molecules of 
length L ranging between 401 and 2080 base-pairs (bp; 
1 bp = 0.34 nm). The DNA is semi-fiexible since L is 
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FIG. 1: Snapshot of TPM Monte-Carlo simulation: The tethered 
DNA molecule is modeled as a coarse-grained chain of A'^ connected 
spheres (various hues) , fixed to the glass surface at one end and to 
the tracked particle at the other end. Here, N = 80, the DNA 
contour length is L = 2080 base pairs (~ 700 nm) and the particle 
radius is B. = 150 nm. 



on the order of the persistence length, £p ~ 147 bp. To 
this end, we combine TPM experiments and dynamical 
Monte Carlo simulations that take into account hydro- 
dynamic effects in the vicinity of the surface. Comparing 
both data sets, we demonstrate that the surface and the 
particle do not affect significantly the polymer dynam- 
ics provided that the particle radius R remains small as 
compared to the DNA contour length. Quantitatively, 
for the polymers considered in this work, this amounts 
to i? < L/6. A radius i? < 20 nm satisfies this con- 
dition for all DNA lengths usually studied by TPM. In 
addition we address rigorously a critical instrumental is- 
sue: the detectors used in TPM experiments always have 
a finite exposure time, ranging from milli-scconds to a 
fraction of second [3 [TOl fTTHSni . We show that exper- 
imental studies must take into account the finiteness of 
this exposure time in order to extract valuable relaxation 
times. Finally, we discuss the intrinsic time resolution of 
TPM experiments when monitoring DNA conformational 
changes and give an illustrative numerical example in the 
case of DNA looping/unlooping. 



formed using a protocol similar to the one described in Ref. [7]. 
A coverslip flow chamber (30 fj,h volume) was incubated with the 
anti-digoxygenin antibody (20 mg/L; Roche) in phosphate-buffered 
saline (PBS) for 20 min at room temperature. After washing, the 
chamber was incubated with casein (1 g/L) in PBS at 4°C for 4 
hours. Experiments with the -R = 150 nm particles (Anti-dig flu- 
orescent particles, Indicia) were performed using chambers whose 
surfaces were derivatized with a mixture of polyethylene glycol and 
biotinylated polyethylene glycol (Sigma, Nanocs) based on the pro- 
tocol of Ref. |21| . The chamber was incubated with neutravidin 
(20 mg/L, Molecular Probes) in PBS for 20 min at room tempera- 
ture. The chamber was subsequently rinsed and incubated with the 
mix of DNA and neutravidin (20 nm, 100 nm) or anti-dig (150 nm) 
coated particles (ratio 1:1, 1 pM) that had been prepared 1 hour 
before in PBS with 0.1 g/L BSA. Due to the lack of small particles 
coated with antiDig, we could not use the PEG passivated cover- 
slips, known to reduce unspecific binding of particle-DNA complex 
on the glass, in all the conditions. 

The tethered particles of 150 nm (resp. 100 nm and 20 nm) were 
visualized by fluorescence video-microscopy, with a magnification 
of 63x4x, on a CCD camera Coolsnap (resp. Cascade, RoperScien- 
tific) at a recording time of 25 frames/s (resp. 74 frames/s and 112 
frame/s), which corresponds to an acquisition period Tac equal to 
40 ms (resp. 13.5 ms and 8.9 ms). In all cases, we could restrict the 
exposure time Tex to 5 ms by the use of an AOTF, characterized by 
an extinction ratio exceeding 7400 (or by entering such a command 
in the acquisition parameters of the Cascade CCD). 

The two-dimensional projection ry of the particle position was 
determined on successive images as the spot centroid using a home 
built image analysis program (Labview) . The accuracy of the posi- 
tion detection (pointing error), calculated as the standard deviation 
of the positions of immobilized particles accumulated during 30 s, 
is equal to 36, 21, and 8 nm for the particles of radius 20, 100, 
and 150 nm, respectively. Trajectories have an average duration 
of 90 s. Along a trajectory, the DNA anchoring point at time t 
is determined by averaging the particle position over an interval 
of duration Tav = 2 s centered at t, and then subtracted from 
ry. This sets (r||) = and subtracts the instrumental drift. Since 
Tav is much larger than the diffusion relaxation time (see below), 
this anchoring point is determined with a good accuracy. The tra- 
jectories exhibiting asymmetry were discarded following 1221 . The 
amplitude of movement of the particle is defined as the variance of 
rii averaged on the trajectory as 



(1) 



A minority of trajectories appeared to have an amplitude of move- 
ment Ary significantly shorter than the majority ones, which yields 

a bimodal distribution of Ar y with a small population of low am- 
plitudes. These trajectories, that might be related to multi-DNA- 
particle complexes, have been discarded. The final number of tra- 
jectories for each condition (L, R, Tex, Tac) ranges from 19 to 60. 
Details are given in the Supplementary Data. 



MATERIALS AND METHODS 
Tethered Particle Motion (TPM) Experiments 

DNA substrates were obtained by PCR amplification from plas- 
mid templates with a 21-digoxigenin-modified forward primer and 
a 21-biotin-labelled reversed primer (Eurogentec) as described in 
Ref. 7 . The DNA substrates DNAR401, DNAR798, DNAR1500, 
DNARL2080 were produced using pAPT72 as a template (posi- 
tions: 1460-1861, 1063-1861, 361-1861, 4625-1861, respectively). 
Their lengths are L = 401, 798, 1500 and 2080 bp respectively. 

Experiments with fiuorescent latex particles of radii 20 and 
100 nm (Fluospheres Neutravidin, Molecular Probes) were per- 



DNA coarse-grained model 

The labeled DNA polymer is modeled as a chain of connected 
small spheres of radius a, whose positions arc denoted by ri{t) 
where i = 0, . . . , N — 1, and a larger final particle of radius R> a, 
of position r{t) (see figures^ and S2). The DNA contour length 
is L = 2a(N — 1). In this work, a ranges between 1.4 and 7 nm, 
and R between (i.e. no particle) and 200 nm. The internal 
structure of the double-stranded DNA is not considered at this 
level of modeling. Denaturation bubbles arc too scarce at room 
temperature to have an effect on the global chain conformation 1231 
I24| . The persistence length value tp = 147 bp is averaged over the 
nucleotide sequence. Torsional degrees of freedom are omitted as 
a first hint into the full problem. The polymer is grafted on a 
surface which sets fq = 0. We use cylindrical coordinates (p, tp, z): 
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r = pBp + ZBz = ry + ze^. Since the polymer motion is limited to 
the upper half-plane, we impose the following "hard wall" boundary 
conditions: Zi > for monomer spheres and z > R — a for the 
particle. We treat a freely rotative joint to the glass coverslip [8], 
by fixing the first sphere center at a height a above the surface (see 
also the Supplementary Data). 

All spheres interact via stretching and bending forces: the po- 
tential U is the discrete version of the extensible worm-like chain 
potential and depends on the sphere positions and the particle 
one r as [/ = [i^ d^. " - 2af ^ (1 - cos 6.)] + 

^(|riv— 1 — r\ — R — 1)^ + ^(1 — cos6jv— i) where Oi is the an- 
gle between neighbouring bonds of sphere i. The first term en- 
sures the polymer connectivity and the second term is the bending 
energy with zero spontaneous curvature. The last two terms are 
dedicated to the particle. The parameters 7 and e are the stretch- 
ing and bending moduli I25| . The persistence length is given by 
£p = e/{k^T). We choose ^a? = ie, which is exact for an isotropic 
elastic cylinder with radius a. Therefore the key parameters in the 
simulation are L/lp, R and A'^. Mutual penetration of monomers 
is prevented by an excluded volume interaction. 

A fully realistic description would be to fix the sphere radius, a, 
equal to the DNA half-width, i.e. about 1 nm. Hence one sphere 
would model roughly 6 bp and 25 spheres would correspond to the 
typical DNA persistence length of 147 bp at physiological tempera- 
ture and salt concentration. Time limitation in the simulations led 
us to concentrate on Af = 50 (or 25) for a DNA of 400 to 2000 bp 
and thus to choose a > 1 nm. As far as equilibrium properties 
are concerned, the DNA statistical mechanics are insensitive to the 
choice of A'^, provided that a remains small as compared to <?p, so 
that a sphere represents a DNA segment that is actually rigid. 



Diffusion near the surface: hydrodynamic effects 

The motion of a spherical particle is slowed down near a flat 
surface due to the no-slip condition for the solvent velocity flow 
at the wall. The induced hydrodynamic interactions cause a vari- 
ation, with the distance to the surface, of the diffusion coefficient 
as compared to its bulk value Dq. This new diffusion tensor can 
be split into a component parallel to the wall, Dy, and a perpen- 
dicular component, D^. For a sphere of radius b (here 6 = a or 
i?), the center of which is at a distance h from the wall, parallel 
and perpendicular diffusion coefficients are derived from Faxen's 
law ^28 . At order 3 in b/h, 

Monte Carlo simulations must be modified to take into account 
these spatially varying diffusion coefficients. First of all, random 
moves 5r are now randomly chosen in an ellipsoid to account for 
anisotropy. In addition, careful attention must be paid to the dis- 
cretization of the equations of motion in this case. A vertical drift 
term (dD^/dh)ez must be added to compensate the variation of 
Dx with h and to restore the detailed balance condition |29) . As 
a consequence, this improvement of dynamics does not affect equi- 
librium properties such as chain statistics. 



Extracting relaxation times from experimental and 
numerical data 



Dynamical Monte Carlo simulations 

Out-of-equilibrium dynamics can be tackled numerically by Dy- 
namical Monte Carlo (DMC) simulations that become equivalent to 
Brownian dynamics when the variation of energy at each time step, 
AC/, satisfies AC/ <IC k^T |26| . We have also performed Brownian 
dynamics simulations of the system which yield very similar re- 
sults on equilibrium and dynamical properties (see Supplementary 
Data). At each Monte Carlo Step (MCStep) of physical duration 
<5i, a bead is chosen uniformly at random among the A' + 1 possi- 
ble ones (monomer spheres and labeling particle). Then a random 
move (5r is attempted for this bead, uniformly in a ball of center 
and radius R^,, thus (6r2) = i-K J^>' r^dr/A-K J^i" r^dr = (3/5)ij2. 
This quantity must be equal to QDoSt, where Dq is the diffusion 
coefficient of the spherical bead, depending on its diameter. In 
practice, for monomer spheres, Ri, = a/5 (unless stated differ- 
ently), where a is their radius. Then St is set subsequently through 
5t = Rl/{10Do), and, for the particle, Rj, = i/IODpt is fixed with 
Dg = Doa/R. Interactions between adjacent beads are treated via 
the interaction potential energy U, whereas interactions between 
non-adjacent beads are of hard core nature, like surface-bead in- 
teractions: whenever a move would lead to the penetration of a 
bead into an other one or the surface, it is rejected. A Monte Carlo 
Sweep (MCS) is a sequence of N + l MCSteps. The physical time is 
incremented of St following each MCS |26| . Typically, a simulation 
lasts between 10^ and 10"'^'^ MCS, which leads to satisfactory error 
bars, estimated by usual techniques [iB\ [see also equation l |12| ]. 
In the simulations, Af = 50 or 25, depending on the slowness of the 
physical process. Since a <g; £p in all cases, we expect that dynam- 
ics do not depend on N either. A simulation snapshot is shown 
in figure^ Finally, note that in our DMC simulations, moves are 
local (one bead at once). This choice is more time-consuming than 
the global Monte Carlo sampling of Refs. [9| |27l. but has the ad- 
vantage to give access to dynamical properties, which is the main 
goal of the present work. 



The relaxation time ry associated with ry can be defined through 
the two-time correlation function averaged over a trajectory 

C{t) = (ri|(s-f t)r||(.)>, ^ (rjj>exp(-t/r||), (3) 

if one assumes without loss of generality that (r||) = 0, or through 
the 2D Mean Square Deviation (MSD): 

MSD(t) = ((r||(s-ft)-r||(s))2), = 2(rj|>-2C(t) 

^ 2(rf>[l-exp(-t/r||)]. (4) 

At short times t <^ ry , one expects to recover the 2D diffusion 
law MSD(t) = 4Dt, with D the apparent particle 2D diffusion 
coefficient. Thus the correlation time for this particle in a 2D trap 
with variance Ar^ is taken to be 1171 

Arjl 

^11 = -■ (5) 

II 2D ^ ' 

For simulated trajectories, the relaxation times are fitted from 
C{t) and MSD(t), using equations | |3|4^ , leading to two relaxation 
times, TQ and tmsd- Below, we report the mean values Tm = 
(tc -|-tivisd)/2, with error bars taken as |t(j — tmsd 1/2. In practice, 
fits are performed on an interval t £ [0, tsup]- Since numerical error 
bars on C{t) and MSD(t) are larger and larger as t grows, the 
smallest possible value of isup must be used. On the other hand, 
for MSD(f), one must have t^up larger than a few tmsd in order to 
fit properly the exponential decay. We have chosen tsup = 4tmsDi 
which is a good compromise between both constraints. The fitting 
procedure, consisting in fitting ttvisd on [0, tsup] and then adjusting 
tsup = 4TMSD! is iterated a few times until tmsd is converged. 
Similarly, is obtained by measuring the slope of ln[C(t)/C(0)] ~ 
—t/rc, on the interval [0,tc], with the same iterative procedure. 
However, in some instances {R < 20 nm), ln[C(t)/C(0)] appears to 
display a short transient, equal to a small fraction of rp, because of 
slow diffusion modes. In this case, linear regressions are performed 
on a suitably chosen interval [tinfiTp], where they appear to be 
very good (correlation coefficients \r\ > 0.9995). 
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For experimental data, ry is fitted from C'{t) as follows. The 
raw correlation function averaged over all available trajectories is 
denoted by Craw(i)- The systematic pointing error is taken into 
account: the detected position, ry is the sum of the actual 
position, ry and the pointing error, Fc, two independent random 
variables. Thus Ar^ = (r?, ) = (r^> + (r^). The sec- 

II, raw \ II, raw' II 

ond contribution is systematically subtracted from measured values 
Ary = Craw{0), using the pointing error values as given above. 
This modified correlation function is denoted by Cm(i). In addi- 
tion, the subtraction of drift induces systematic anti-correlations 
at short times leading to the following fitting form 



Cm(f) 
Cm{0) 



H-2 



Tm 



-t/T 



(6) 



with a single fitting parameter, Tm. Examples are given in figure SI. 
The prescription is the same as above: the fitting interval [0, tsup] 



is chosen so that ta 



4Tin. Error bars on Tm correspond to the 



standard deviation of the measurements on individual trajectories. 



Correction of detector time-averaging effects 

Finally, one has to correct time-averaging effects in experimental 
results. In Refs. (9] 1101 117lt20| , the time-averaging (or blurring) 
effect due to the finite exposure time of detectors in single molecule 
(or particle) tracking experiments was investigated. When tracked 
molecules or particles diffuse in confined regions, diffusion constants 
can be significantly under-estimated, as well as sizes of confining 
domains. This effect was quantified by exact analytical arguments, 
in the contexts of diffusion in membrane domains 1171 and optical 
traps 1181 . leading to the same correction to Ary. In the present 
work, the situation is similar because the motion of the particle is 
restricted by the DNA tether, thus all the analytical derivation of 
Refs. [171 118| remains valid: the real relaxation time ry , domain 
size (or root-mean-square excursion from the attachment point) 
Ary and diffusion coefficient D can be related to their measured 

counterparts, r^, Ary ^ = \/C\-n{0) and D^. Whenever > 

2 Tox/3 [la 

Ty ~ Tm - Tex/S, (7) 

where Tex is the detector exposure time, and 



Aril 



: Ar,, 



1 -1/2 



(8) 



If Tni < 2 Tcx/3, no correction can be applied |17l and one ought 
to switch to a faster acquisition device. 



RESULTS AND DISCUSSION 
Equilibrium distributions of particle positions 

Before considering the dynamics of the particle-DNA 
complex, we first briefly analyze the particle positions at 
equilibrium through the histograms of projected particle 
excursions, p = jryl, i.e. the probability distributions 
p{p), and through the standard deviations Ary = (r|)^/^. 
We compare the latter to the fitting formula suggested 
by Towles et al. \W\ , in order to test both our simulation 
and experimental results. 

We simulate DNA molecules of lengths L = 
401,798,1500 and 2080 bp, corresponding to the semi- 
flexible case (2 < L/£p < 14) for which no analytical 
expression for p{p) is known. The theoretical study of 
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FIG. 2: Comparison between experimental (solid symbols) and 
numerical (open symbols) standard deviations Ary = (rl)-*^/^. Ex- 
perimental values have been corrected for systematic errors. Error 
bars are 95 % confidence intervals. Red and red dashed lines are 
predictions from Ref. [9] for = 147 and 128 bp, respectively. 



the limiting cases, rigid rod {L ^ £p) and flexible chain 
(L ^ £p), is reported in the Supplementary Data. In ad- 
dition to the experimental particle radii R = 20, 100, and 
f 50 nm, we also examined numerically the radii R = 40, 
80, and 200 nm, as well as i? = (no particle). The 
numerical distributions p{p) and their deviation from a 
Gaussian are discussed in the Supplementary Data (fig- 
ures S4-S6). 

The two distributions p{p) extracted from both TPM 
experiments and simulations are compared in the Supple- 
mentary Data (see figures S5, S6). The agreement is sys- 
tematically good, even though some small discrepancies 
appear. In order to quantify them, we plot in figure[2]the 
values of Ary as a function of L. Interpolating functions 
evaluated from Monte Carlo simulations have been pro- 
posed in Ref. 9 that give Ary in function of L, R and 
£p. We have checked that our numerical results are in 
agreement with these interpolations when R > 100 nm, 
even though they were calibrated for R > 190 nm. In 
figure [2j one observes that experimental Ary are shorter 
than simulated ones. Nelson et al. [9] have already no- 
ticed that assuming the persistence length £p = 147 bp 
(50 nm) in data analysis leads to larger Ary than ob- 
served. To circumvent this issue, they pointed out that 
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setting £p = 128 bp (43.4 nm) yields better agreement, 
attributing this low value to the buffer. We also observe 
that setting £p = 128 bp leads to values of Ary closer to 
our experimental observations for R — 150 nm. 

However, for i? = 20 nm, the interpolating functions 
proposed in Ref . |9] are less reliable (see figure [2]) . In 
particular, decreasing the persistence length value to 
£p = 128 bp is no longer adequate and goes into the 
wrong direction (147 bp is a better choice). 

Note that the actual value of £p is certainly not the 
unique source of discrepancy between theory and exper- 
iments. For example, the interaction between the func- 
tionalized surface on the one hand, and the particle (and 
DNA) on the other hand, is more complex than a simple 
excluded volume one. Therefore inferring a precise value 
of £p from such experiments is a challenging task [TTl |3D] . 
Given the approximations used in the physical modeling, 
we conclude that the agreement with experiments (dis- 
crepancies are less than 15%) is satisfactory. 



Relaxation times of the particle-DNA complex 

In this section, we compare the relaxation times ry 
extracted from both experimental and numerical trajec- 
tories (see Methods) . We first present theoretical consid- 
erations on dynamical properties of a DNA fluctuating 
freely in solution and then give basic insights into the 
influence of the attached particle. We finally address the 
issue of the whole particle-DNA-surface system, both ex- 
perimentally and numerically. 

Theoretical aspects on relaxation times and diffusion 
coefficients - For a free chain of polymerization index N 
and radius of gyration Rq, the correlation time is |31j 



D. 



chain 



Do ■ 



(9) 



In the second equality, we have used the Rouse model 
ignoring hydrodynamic interactions and valid for flexible 
chains {L/£p ^ 1), which sets that Dchain = Dq/N where 
-Do = ksT/ {GTTrja) is the Stokes' formula for the diffusion 
coefficient of a monomer sphere of radius a in a liquid of 
viscosity rj. For a Gaussian chain, Rq cx N and we get 
Ty ^ N'^/Dq at large N. For a rigid rod, i?g cx N'^ and 
ry ^ N^Do. 

Note that at the level of modeling chosen in the present 
work where the real polymer is modeled by N connected 
monomer spheres, dynamics (and thus ry) is not affected 
by the choice of N as soon as a ^ £p. Indeed, Rg being 
insensitive to the choice of a (or iV) , the only dependence 
in N comes from the ratio N/Dq in equation ([9]). But 
with the definition of Dq, N/Dq cx Na cx L, indepen- 
dently of a. 

How is this chain relaxation time modified when a par- 
ticle is grafted at one end? The time evolution of the 
DNA polymer and the particle is governed by the over- 



damped Langevin equation [3T] 

CMt) - -Vr,C/(ro,...,rw-i,rpart)+e.(t), (10) 

where are the position of spheres {i = 0, . . . , iV — 
1) and the particle {i = part) (we neglect hydrody- 
namic interactions as a first hint into the full prob- 
lem). This equation relates the linear response of each 
object i to forces applied to it: on the one hand, the 
derivative of the potential U (sum of pairwise poten- 
tials) between the N objects; on the other hand, the 
stochastic forces ^i{t) mimic the action of a thermal 
heat bath and obeys the fluctuation-dissipation relation: 
{ii.it) ■ ^j(^')> = 6kBTQS{t-t')5ij, where the friction co- 
efficient Ci {— Co for monomer spheres and (^oR/a for the 
particle) is related to the diffusion coefficient Di through 
D, = kBT/Q. 

We relate the diffusion coefficient of the particle-DNA 
complex, Dc, to those of the particle alone, -Dpart, and 
the DNA polymer alone, Dchain- To get rid of the inter- 
action forces, we consider the barycenter of the particle- 
DNA complex Tc = (Cchainfchain + Cpartrpart)/Cc whcrC 

Tchain = DNA Center of mass. Then 

a linear combination of Langevin equations ( 10 1 yields 
Ccrc(0 = ic{t), where |c 



schain 



spart ■ 



The correla- 



tion function (^c(0 ■ lc(i')) = GkBT^J^t - t') sets the 
diffusion coefficient of the particle-DNA complex, Dc- 



d: 



D 



part 



D 



chain ' 



(11) 



We emphasize that Dc should not be confused with the 
measured effective particle diffusion coefficient, D, be- 
cause Tc 7^ fpart- The only case when they are equal is 
for -Dpart ^ Dchahi- By contrast, the relaxation time ry, 
obtained by tracking the particle only, is a feature of the 
dynamics of the whole complex. This is the reason why 
we focus on this observable in the following. 

Simulation results without hydrodynamic interac- 
tions - We first display in figure [3] the relaxation times 
fitted from numerical trajectories simulated without hy- 
drodynamic interactions with the surface. They are de- 
noted by Ty™. We checked that finite- effects appear 
to be negligible for a ^ £p (data not shown). We also 
explored the reference case i? = 0, where no particle is 
attached to DNA. We have seen above that the Rouse 
model predicts ry ^ at large L. Figure [s] shows that 
this regime is also a good approximation even for the fi- 
nite sizes considered here. We also plot in figure |3] the 
relaxation times in the case where i? = in the absence 
of hard wall at z = 0. One can see that the wall does not 
hinder DNA dynamics significantly. 

When R > 0, the particle slows down dynamics be- 
cause viscous drag increases with the particle size. How- 
ever our simulations predict that DNAs are not signif- 
icantly affected by a small particle {R = 20 nm). A 
medium particle, i? = 40 nm (resp. 80 nm) does not 
perturb DNAs of lengths L > 800 bp (resp. 1500 bp). In 
the other cases, increasing R from to 150 nm at fixed 



slope 2 

1000 

L(bp) 

FIG. 3: Numerical relaxation time Ty'™ (solid lines: without wall 
hydrodynamic friction (no HI); N = 50, except for L = 401 bp, 
R = 150 nm and L < 798 bp, _R = 80 and 200 nm, where TV = 25) 
and ''"ii™ (dashed lines: with wall hydrodynamic friction; N = 25 
for L = 401 and 798 bp, AT = 50 for L = 1500 and 2080 bp) 
versus the polymer length L for various particle sizes R, in log-log 
coordinates. In the case R = 0, Ty is also plotted when there is no 
hard wall. The dotted line of slope 2 shows the expected scaling 
for long DNAs (Rouse model for a random walk). 
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FIG. 4: Experimental (''"|^''^i solid symbols) and numerical (r^'™, 
using z-dependent diffusion coefficients, open symbols) relaxation 
times for different DNA lengths, L, and particle radius, R, in linear- 
log coordinates. Displayed error bars are explained in the text for 
numerical data and are 95 % confidence intervals for experimental 
data. Other error bars are smaller than symbol sizes. Inset: Ratio 
T|°''''/t|™ versus L/R with the same symbols as above, in log-log 
coordinates. The dashed line shows the best linear regression, with 
slope -0.355 ~ -1/3. 



L monotonically increases Ty™, by a factor ranging from 
2 {L = 2080 bp) to 10 (L = 401 bp). 

In the Supplementary Data (see figure S8), we com- 
pare our numerical values of figure [3] without wall hydro- 
dynamic friction to simple calculations where the DNA 
molecule is modeled by an ideal spring. In this case, the 
relaxation time would be equal to the ratio of the parti- 
cle (or particle-DNA complex) friction coefhcient to the 
DNA spring constant [TT]. This approach appears to be 
insufficient because the interplay between the polymer, 
the bead and the wall is more complex than this simple 
image. 

Simulation results with hydrodynamic interactions with 
the surface - Hydrodynamic interactions between the 
particle-DNA complex and the surface are implemented 
in the numerical code using equation (2]). Figure [s] shows 
that the diffusion is slowed down by hydrodynamic cor- 
rections for large particles {R > 80 nm) but is weakly 
affected for small ones {R < 40 nm). Indeed, small 
particles are far away from the surface, in that sense 
that R z and b/h = R/{z + a) is generically small 
in equation Only the first DNA spheres, for which 
b/h = a/{zi -\- a) is of order 1, are slowed down by the 
proximity of the surface but this is not sufficient to affect 
the whole dynamics. This is related to the observation 
that the exact modeling of the DNA-surface joint is not 
a relevant issue (see Supplementary Data). By contrast, 
when R is large, R/ {z + a) is close to 1 for the particle 
itself, and hydrodynamic corrections play a significant 
role. In the following, we use these numerical results 
with hydrodynamic interactions. 



Comparison between experimental and numerical re- 
laxation times - Figure |4] displays our experimental [cor- 
rected from detector averaging effects using equation ([7|] 
and numerical values of the relaxation times ry extracted 
from functions C{t) or MSD(i) (see Materials and Meth- 
ods). Experimental and numerical values are found in 
good agreement, with ratios of experimental to numeri- 
cal values varying from 0.5 to 2. These ratios are plotted 
in the Inset of figure [4j They appear to be well correlated 
to the ratios L/R and to follow the approximate power- 
law: T-|P/r||'" cx [R/Lf/^. This suggests that the ob- 
served discrepancies are not due to statistical errors (as 
supported by the small error-bars), but have a physical 
origin. When R ^ L, t^^^'^ < i.e. experiments are 

faster than simulations. If i? > L/5, then t^'''^ > 
Indeed, if one neglects hydrodynamics, the measured dif- 
fusion coefficient of the particle-DNA complex, Dc, is re- 
lated to the particle one, Dpart, and to the polymer one. 



i^chain, through equation (|11|: ^ = Dp^^^. + f^hlin- 
Thus the particle does not slow down the complex pro- 
vided that Z3part > i'chain- Now impart ^ K/R with 

K = kBT/iQmf), and Dchain = K/{Na) = 2K/L in the 
Rouse approximation. If 2R <^ L, then Dpait ^ -Dchain 
and the DNA dominates the particle-DNA dynamics, as 
if the particle were absent. But in this case, hydro- 
dynamic interactions between DNA segments should be 
taken into account. Since such interactions are known 
to accelerate polymer dynamics [25l[3T], one expects ry 
to be shorter than in simulations, as indeed observed in 
experiments. Conversely, if 2R ^ L, then the particle 
dominates the particle-DNA dynamics, as if the role of 
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the DNA polymer was limited to its spring properties. 
We have simplified the complex hydrodynamics of this 
particle in the vicinity of a wall through the expansions 
in equation ([2|. But these expansions are valid in the 
small b/h limit only and we cannot expect them to be 
correct when b/h ~ 1, which is the case when 2R ^ L. 
In Ref. [23] Section 7-4, it is proved that equation ([2]) 
overestimates D_l and Dy in this case. Consequently, 
T|| is underestimated as observed. Describing correctly 
the particle dynamics would require a full integration 
of all hydrodynamic interactions, between the monomer 
spheres, the particle and the surface. This is beyond the 
scope of this work. 



Choice of TPM parameters to infer DNA dynamics 

The preceding analysis leads us to suggest the following 
recommendations to extract the most valuable informa- 
tion on DNA dynamics from TPM experiments. 

Small particles of radius i? < 20 nm do not slow down 
DNAs as short as 400 bp - Comparing simulation data 
for i? = and R > (figure [s]) shows that the pres- 
ence of the particle does not increase the relaxation times 
for small particles or long DNAs. Indeed, provided that 
R ^ L/2, dynamics are dominated by the polymer. Ac- 
cordingly, figure [sl shows that if i? < L/6 then ry is less 
than twice the relaxation time of DNA without particle, 
which demonstrates that slowing-down is weak or mod- 
erate in this case. In particular, slowing down is very 
weak when i? = 20 nm for L = 401 to 2080 bp. 

They are three contributions to hydrodynamics: i) 
slowing down of the chain and the particle, due to the 
no-slip condition at the surface, ii) hydrodynamic in- 
teractions between the chain and the particle, and iii) 
interactions between monomer spheres. For creeping 
flows, these three contributions are additive [55]. We 
have shown above that, for small R, (i) affects only very 
weakly the relaxation times. Contributions (ii) and (iii), 
which are not taken into account in simulations, moder- 
ately accelerate the dynamics of the particle-DNA com- 
plex, as corroborated by experimental measurements (di- 
amonds in figure [4]) . By linearity, we thus expect that 
neither the wall nor the small particle slow down DNA 
dynamics as compared to a DNA chain fluctuating freely 
in solution. 

Taking into account exposure times and acquisition pe- 
riods - In the Materials and Methods section, we have 
explained how the flniteness of the exposure time, Tqx, 
modifles the root-mean-square excursion of the particle, 
Ar|| , the particle-DNA relaxation time, ry , and the ef- 
fective particle diffusion coefficient, D (when Tox < ^ry). 
We have prescribed how to recover their real values rig- 
orously, contrary to the rough approximation of equa- 
tion ([s]) proposed by Ref. [TUI. Indeed, the formula given 
in this anterior work largely overestimates the averaging 



L(bp) Tex(ms) 


Tm (ms) 


Ary,m(nm) 


ry (ms) 


Ary (nm) 


2080 40 


69.5 


223 


56.1 


250 


5 


59.1 


258 


57.4 


261 


401 40 


32.8 


91 


19.5 


122 


5 


21.7 


118 


20.0 


123 



TABLE I: Relaxation times and movement amplitudes before (rm 
and Ary j„) and after correction (rn and Ary) of detector time 
averaging effects, using equations ( |7|8| , for two DNA lengths L and 
two exposure times Tex. Here the particle radius is R = 150 nm. 



effect, thus requiring the introduction of a phenomeno- 
logical time scale t* . 

To this respect, we have analyzed two particle-DNA 
complexes {R = 150 nm and L = 2080 and 401 bp) with 
two different exposure times, Tox = 5 and 40 ms and the 
same acquisition period Tac = 40 ms in both cases. The 
measured values before and after applying equations ( 7p ) 
are given in Table |T] We observe that the corrected values 
are in excellent agreement for both DNA lengths. These 
values validate our method to correct detector averaging 
effects the experimental data. 

The values of ry given in figure |4] will help experimen- 
talists to anticipate the appropriate required values of the 
exposure time Tgx in future works. Owing to the com- 
plexity of the problem, no simple formula can be given 
that expresses ry as a function of R and L. In practice, a 
linear interpolation of the values of figure [4] will provide 
a good estimate of ry for other values of R and L, in the 
range studied. 

However, reducing Tox may remain insufficient when 
one is interested in measuring ry . The acquisition period 
Tac being imposed by the camera, the only points avail- 
able to fit MSD(i) (X l-exp(-i/ry)orC(i) (xexp(-i/Ty) 
plots are measured at discrete values t — 0, Tac, 2Tac and 
so forth. So if Tac is large as compared to ry, MSD(Tac) 
and C(Tac) have already reached their asymptotic values, 
respectively 2Ar| and 0, making any fitting procedure 
hopeless. If one is interested in dynamical properties, 
the true technical limit obviously remains the camera ac- 
quisition period, which must be of the same order of mag- 
nitude as r,„ or smaller, typically Tac ^ 2rni (as detailed 
in the Supplementary Data). 



Measuring variations of (r|) in real time - Beyond the 
measurement of equilibrium statistics or relaxation rates, 
TPM experiments can also be used to monitor DNA con- 
formational changes. To what extent can those confor- 
mational changes be monitored in real time? In TPM, 
such a change is characterized by a variation of the move- 
ment amphtude, Ar| = (r|), due to a given molecular 
event such as protein binding or DNA looping, leading 
to a variation of the apparent DNA length. But mea- 
suring the mean value (r^) supposes to average r| on 
a sufficiently long time interval, of duration again de- 
noted by Tav [I]. By definition of the relaxation time 
ry , the average will be accurate provided that Tav ^ '''y • 
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Once ry is known, what averaging time Tav should be 
chosen in practice? In other words, what is the mini- 
mum time scale of conformational changes that TPM can 
reach? Formally speaking, (r|) ~ J^'j'^^'"' r|(i) dt/Tav 
The relaxation time of r|, T|| 2, is again defined by 
(rj|(s + t)r|(s))5 - (r|>2 ss cr^^ exp(-t/T||^2) (we use the 
notation 0^.2 for the r.m.s. of rn, to avoid confusion with 
Ar|, the variance of r||). Then the statistical error on 
the measurement of (r|) (68% confidence interval), de- 
noted by Err(ry), can be quantified as follows 26 : if 
T\\,2 ^ T'ac, then Err(r^) = i/Sry^aT^crrij • This rela- 
tion takes into account statistical correlations between 
successive frames, measured by T|| 2- If th 2 ^ Tic, then 
Err(r|) = ^JT^^^T^o ^1 ^ because all frames are statisti- 
cally independent. 

Now, what is the relationship between ry 2 and ry? 
Since = a;^ + j/^, T|| 2 = = Tj,2, the (identical) re- 
laxation times of the obscrvables and . \i x and y 
were Gaussian random variables, then t^^i — Tx/2 = Ty/2. 
But X is not exactly Gaussian. In practice, we have mea- 
sured numerically in our case that ry 2 — ■'"||/c with c 
ranging from 2.5 to 3.5, depending on L and R. Thus 
'''11,2 ^ ''"II /2- Furthermore, we have also measured numer- 
ically that (7^2 < (ry) for the values of L and R studied 
here. Were ry be exactly Gaussian, a simple calculation 
shows that (7^2 = (r|) (because r|| is two-dimensional). 
Thus in the regime where rn 2 ^ Tic, 



is to show how the minimum detectable lifetime can be 
determined a priori from the knowledge of rn . 



A = 




(12) 



The equality defines the limiting value T*^ = t\\/\^. For 
example, if Tac = 10 ms and T|| w 30 ms, Tav > T^av — 
100 ry 3 s yields a relative statistical error A < 10%. 
If the conformational changes lead to variations of (r|) 
larger than 10% (the statistical noise), then a threshold- 
ing method [4l El [131 E] can detect them in real time 
(provided that their time scale is higher than 3 s). Using 
shorter DNAs, smaller particles, and a faster acquisition 
device improves the time resolution. Concrete illustrative 
examples are given in figures S9, SIO and thereafter. 

Numerical example 0/ conformational changes detec- 
tion: DNA looping /unlooping - To further illustrate 
TPM capabilities in terms of detection of conformational 
changes, we propose the following example. It is re- 
lated to TPM experiments where DNA loops are pro- 
moted by proteins binding to two specific DNA loci, and 
DNA conformation alternates between looped and un- 
looped states [H [SI US]. In the looped state, the effec- 
tive DNA contour length is shorter, resulting in a smaller 
Ar| . Switching between looped and unlooped states can 
be detected in real time provided that their lifetime is 
long enough. One of the purposes of the present example 
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FIG. 5: Simulated amplitude of movement Ar|(t) for L = 
798 bp, i? = 20 nm and three averaging times Tkv = 3 (yel- 
low), 30 (red) and 300 ms (black). The vertical lines indicate 
the looping and unlooping events in the simulation. The hor- 
izontal solid lines show the expected values of Ary in the 
looped (bottom) and unlooped (top) states. The dashed line 
indicates the threshold separating these two states for detec- 
tion purposes. The bottom plot is a zoom of the top one on 
a time interval where conformational changes are very close. 



We focus on L = 798 bp and i? = 20 nm with TV = 25 
beads (denoted by Bq, . . . , i?24)- Hydrodynamic inter- 
actions with the wall are taken into account. The nu- 
merical relaxation time is ry = 2.98 ~ 3 ms (see fig- 
ure l3|, and Ary = 124 nm. The two loci where pro- 
tein oinding occurs are situated on beads and Big. 
Thus the looped DNA is roughly equivalent to a 400 bp 
long DNA, with Ary = 81 nm. In the simulation, 
looping is promoted by an attractive quadratic poten- 
tial between beads and B\^: V — 2fcBT"(ri8 — rg)^. 
When we turn on (resp. off) the potential, the DNA 
molecule switches into the looped (unlooped) state. In 
figure jsj Ar^ is plotted in function of time for differ- 
ent averaging- interval lengths Tav To set the threshold 
separating looped and unlooped states, we choose the 
condition (H-A)Ar| j^^p^^ = (1 - A)Ar|^^i^^p^j, leading 
to A = 0.4. This means that a maximum relative error 
of 40% is allowed if one wants to detect unambiguously 
looped and unlooped states. Using equation (12), this 



leads to the limiting value T*, 



19 ms (pro- 



vided that Tac ^ T\\.2 ~ 1 ms). As illustrated in figure [5j 
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choosing Tav < T*^ leads to very numerous false detec- 
tions of looping/unlooping. On the contrary, increasing 
Tav above the threshold T*^ makes false detections un- 
likely. However, as is shown in figure [5] increasing Tav 
increases the error on the determination of the moment 
when the transition between looped and unlooped states 
occurs. Hence the looped/unlooped state duration mea- 
sure is less precise (see the zoom) . This precision is on the 
order of magnitude of Tav. Furthermore, because of aver- 
aging, only durations larger than Tav can be accurately 
detected. For instance, in figure [s] (zoom), the looped 
state between t = 5.72 s and t = 5.86 s is not detected 
because its duration is shorter than Tav for Tav = 300 ms. 
Therefore Tav must be large enough to distinguish be- 
tween looped and unlooped states, but it must be as short 



as possible to detect the shortest lifetimes. Equation ( 12 ) 



provides the optimal choice. Of course, the shortest T|| , 
and therefore particles of radius R < T/6, will provide 
the best detection capabilities. Another illustrative ex- 
ample is given in figure S8 in the case of DNA bending 
induced, for instance, by protein binding. 



CONCLUSION AND OUTLOOK 

This work focuses on the influence on DNA dynam- 
ics of the TPM setup geometry, namely the attachment 
of the DNA tether ends to the surface and the labeling 
particle. To what extent DNA characteristic times can 
be inferred from TPM observations? The answer to this 
complex issue cannot be summarized in simple formulas 
because many parameters are involved and interrelated, 
in particular: the particle size R, the DNA contour length 
L, the DNA persistence length £p, the exposure time Tex, 
the acquisition period Tac, and the averaging time Tav. 
Here, we used numerical simulations, the results of which 
were compared to experimental measurements. A good 
quantitative agreement (see figures [2] and [4]) allows us 
to draw the following two major conclusions, summariz- 
ing the survey of experimental constraints on dynamical 
TPM examined in the last section of our Results and 



Discussion. 

1) For DNA lengths L usually used in TPM experiments, 
i.e., from hundreds to thousands base pairs, the DNA 
dynamics is weakly perturbed by particles of radius R < 
L/Q. As long as L is larger than 400 bp, the influence 
of the particle is negligible if i? < 20 nm. Indeed, in 
this case, hydrodynamic friction induced by the particle 
and the surface are weak and the DNA characteristic 
relaxation time in the plane parallel to the surface, T|| , is 
almost identical to the case where the particle and the 
surface are absent. 

2) In the TPM analysis of DNA conformational changes, 
the use of particles of radius smaller than 20 nm allows to 
detect molecular events with a time resolution as short as 
20 ms in favorable cases. We believe that such a substan- 
tial improvement will provide a more detailed view on 
DNA-protein interaction processes involving DNA bend- 
ing [H [15] , looping [H [6l [13] or migration of a Holliday 
junction [5] in the future. We have shown in the simu- 
lation of looping/unlooping events that the knowledge of 
the relaxation time T|| is a prerequisite in order to extract 
the most relevant information from TPM experiments ex- 
ploring DNA conformational changes: if ry is not known, 
it is not possible to decide whether a detected variation 
of Ar| is a real conformational change or a spurious sta- 
tistical fluctuation. 

On the theoretical side, the scaling law t™'' cx 

(i?/T)^/^Ty™ in figure jlj provides a clue that hydrody- 
namic interactions play a crucial role. Brownian dy- 
namics simulations (see Supplementary Data) will be the 
next step in order to elucidate this issue, at least for 
the fastest complexes (short chains and small beads). 
They ought to include hydrodynamic interactions be- 
tween beads [25l |32] , usually taken into account through 
the Rotne-Prager tensor, and hydrodynamic images in- 
duced by the presence of the particle and the surface to 
ensure the no-slip condition on these surfaces |25| . 
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Tethered Particle Motion (TPM) Experiments 



The number of experimental particle trajectories fol- 
lowed in TPM experiments for various DNA lengths L, 
exposure times Tox and particle radii R are given in the 
following table (distances in nm, times in ms): 



DNA length 
L (bp) 


Tex = 40 
R = 150 


Tex — 5 


i? = 20 


R = 100 


R = 150 


401 


46 






32 


798 




34 


24 


32 


1500 




38 


19 


40 


2080 


33 


60 


41 


29 



An additional set of 35 trajectories were followed for 
L = 2080 bp, R = 150 nm, Tex = 5 ms and T^c = 12.8 ms. 



-0.2 I ' < ' < ' ' ' ' ' ' ' ' 

0.02 0.04 0.06 0.08 0.1 0.12 

t(s) 




Fitting procedure 
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The fitting procedure to extract the relaxation time T|| 
from experimental data is presented in the Material and 
Methods section of the paper. The fitting form is given in 
equation (6), which we first demonstrate. The modified 
correlation function C,n(i) is obtained after subtraction 
of noise, that we thus ignore here, and of instrumental 
drift. To subtract instrumental drift, the measured po- 
sition r|[ is first averaged on a time interval of duration 
Tav centered at t: 



Ta. 



av Jt-Tav/2 



r|| {i)di, 



(SI) 



and then subtracted from rn (t). Thus 



CnM = ([r|| {s + t)- r|| {s + t)] ■ [r|| (s) - r,, (s)]),. (S2) 

We assume that r,„ <C Tav and t ^ Tav , which is satisfied 
in practice in this work. We set b — 2rm/rav ^ 1. The 
measured position rn (t) is the sum of the real position, 



denoted by r[°^(t), and of the drift d{t): r\\{t) ^ rf|''^(t)-|- 
d{t). The drift is a slowly varying function of time, so 
that one also has r|| (t) = r["^(<) + d{t). Thus 

Crait) = {[r\"\s + t)- f["'(s + t)] ■ [r\°\s) - f\"\s)]) 



(0), 



(0), 



-Ml 



,(0) 



(0)/ 



FIG. SI: Two fits of experimental correlation functions 
Cm(t)/Cin{0) according to equation (6). Stars are experimental 
points and red curves are fits. Left: L = 1500 bp, R = 100 nm, 
Tac = 13.5 ms and Tex = 5 ms. The fit gives Tm = 25.6 ms. Right: 
L = 2080 bp, R = 150 nm, Tac = 40 ms and Tex = 5 ms. The fit 
gives Tm = 59.0 ms. 



If, without loss of generality, we assume that = 
on each trajectory, then the first average in equa- 
tion (S3) is equal to ([r[°']2)e"*/^"' . Using equation (SI I, 



^r|| '(^sj.r||-'(,ij; = {[rf\-)e-'-^-'" and Tav > t,Tm, one 
gets that the three last averages are equal to fe([r[|'^-']^). 



Therefore C^it) = ([r[°^]2)(, 

J0)l2 



-t/r„. _ ^) and Cm(0) = 
([r^i''^]2)(l - b). It follows that Cm(i)/Cm(0) = (1 + 



fo)e^*/'^" — 6 at first order in &, as claimed above. Fig. SI 
provides two generic examples of such fits of modified 
experimental correlation functions Cm(i)- These figures 
illustrate the quality of the fits. 

An alternative way to circumvent systematic errors on 
Craw(O) consists of excluding Ciaw(O) in the fitting pro- 
cedure and using the fitting formula 



Craw(i) - Arl 



Ta, 



- 2- 



Ta, 



(S4) 
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with two free parameters, Ar^^ p and t^. For all data 
sets, the so-obtained values of r„i are equal, within a 
few percents, to the ones coming from equation (6) of 
the main text, thus confirming the validity of our fitting 
procedure. 



Brownian Dynamics simulations 

For the numerical iterations, the evolution of each 
sphere position f; = r^/a is governed by an iterative 
Langevin equation with a discrete time step St. In terms 
of the discrete time variable n — t/St, 

r,in+l) = t,{n)-DoV,V{{H}) + \ftl)olin), (S5) 

where U = U/ksT and the rescaled random displace- 
ment ^j(n) has variance unity • ^,j{m)) ~ 3(5.y(5„,„ 
[see also equation (10) and below]. The rescaled bare dif- 
fusion coefficient Dq = DoSt/a^ is the diffusion constant 
in an unbounded space in units of the particle radius a 
and time step St: a sphere typically takes 1/Dq itera- 
tions to diffuse on a distance a. For sufficient numerical 
accuracy we choose time steps in the range Dq = 10"'^- 
10~^ to avoid unphysical large displacements and forces. 
Output values are calculated every lO'^-lO* steps, total 
simulation times are on the order of 10^-10^ steps, giving 
reasonable error bars. Finally, since the polymer motion 
is limited to the upper half-plane z > 0, we use the follow- 
ing reflection boundary condition: if a sphere intersects 
the substrate, its height Zi is replaced by z[°**: 



We have observed that, up to error bars, statistical and 
dynamical properties are similar. In terms of acceptance 
ratio T, we have measured that our Monte Carlo param- 
eter choices lead to 68% < T < 88% for all the situations 
studied in this paper, and T > 80% for simulations with 
hydrodynamic interaction with the wall (those compared 
to experiments). To this end, Rt (defined in the Mate- 
rials and Methods section) was set to Rb — a/5 for all 
data sets, with the few following exceptions: Rb = a/7.5 
for L = 798 bp and R = 0; Rb = a/10 for L = 401 bp 
and i? or 20 nm. 

Thus the criterion AU <ti ksT is not satisfied stricto 
sensu, but the fact that both types of simulations give 
comparable results even at the dynamical level proves 
that this value of T is sufficient in our case. The choice of 
excluded volume interaction in the DMC code (whenever 
a move would lead to the penetration of a bead into an 
other one or the substrate, it is rejected) saves computa- 
tional time as compared to the calculation of truncated 
Lennard- Jones potentials used in Brownian Dynamics. 
Since DMC is favorable in terms of computational time 
(about 30 times faster in our case because it saves the 
computation of Lennard- Jones potentials and allows for 
a larger time step), we restrict our numerical work in 
the paper to DMC simulations. Note however that these 
Monte Carlo simulations remain time-consuming: each 
simulation for a set of parameters typically requires sev- 
eral days on a standard processor. 



Exact equilibrium distributions in the rigid and 
flexible limits 



Zi < a 
zjY < R — a 



,rcfi 



2a - z., for i<N (S6) 
2(i? - a) - zjv (S7) 



The excluded volume interaction is modeled by a 
repulsive, truncated Lennard- Jones potential V^j = 
E,<,[(V|f» - f,l)'' - 2(V|f. - + 1] valid for sepa- 
ration jf^ — f j| < 2 and h — R/a + 1 for j — N and 6 = 2 
otherwise. 



The particle center distribution G(ro, r; N) is the poly- 
mer propagator from Tq = to the center of the particle 
r, in N steps (in the z > half-space). The experimental 
observable is not this full distribution but the marginal 
distribution of p 



pip) 27rp 



G{0,r;N) dz. 



(S8) 



R-a 



Comparison between Brownian Dynamics and 
dynamical Monte Carlo simulations 

We have simulated our particle-DNA complex for dif- 
ferent choices of L and R, using both Brownian Dynamics 
and Monte Carlo simulations. 

The relevant test of the dynamical Monte Carlo (DMC) 
algorithm reliability is the acceptance ratio T, i.e. the 
number of accepted Monte Carlo moves as compared 
to the number of attempts, which must be close to 
1. Indeed, in Monte Carlo algorithms, a move is al- 
ways accepted if AU < and is accepted with prob- 
ability expi-AU/ksT) if AU > 0. Thus T ~ 1/2 
(exp(-Ai7/fcBr))Ac/>o/2 ^ 1 - {AU)Au>o/{2kBT) and 
the condition AU ^ fceT is equivalent to T close to 1. 



where Zmax (p) = \/ {Ro + L)^ — is th e ma ximum value 
for z at fixed p and Rq = R + a (Fig. [S2|. This prob- 
ability distribution can be computed analytically in the 
two following limits: rigid rod {L/£p — 0) and phan- 
tom Gaussian chain {L/£p — c»). In the rigid rod limit, 
the propagator is simply fixed by |r| i?o -|- i, in other 
words G'RR(0,r;7V) = A5[y^p'^ + z^ - {Rq + L)] where 
A = [27r(i?o -I- L)(2a + L)]^^ is the normalization con- 



stant. Using equation (S8|, we find 



Prr(p) = 



P 



{2a + L)^{Ro + L)^-p^ 



(S9) 



for < p < p,nax = V(^o + i)2 - (i? - a)2 (FigJS2|. In 
Fig. S3i, we display Prr(u) where u = p/pmax- This dis- 



tribution is linear for small p. When i? ^ L, it sharply 
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FIG. S2: TPM geometry. The particle has radius 7?. The 
polymer, of curvilinear length L, is modeled by a sequence of 
N spheres of radius a. The polar coordinates of the particle 
in the (xy) plane are denoted by (p, (p), pmax is the maximum 
extension and the particle height is denoted by z. The whole 
chain is confined in the upper half-space. For a given p, z < 
Zmax(p). Only p and (p (or x and y) are accessible by TPM. 



increases for p — >■ pmax because of projection effects. 
When R ^ L, this distribution becomes prr{u) ~ 2m, 
owing to Prr(u) = Prr(/o) dp/du = p^ax Prr(p)- 



In the flexible or phantom Gaussian regime, the end-to- 
end distance probability distribution for the chain alone 
is given by G'chain(0, rAr_i; iV- 1) cx exp{-3r%_^/[2{N - 
l)(2a)^]}, where rjv-i is the 3D position of the last 
sphere. The distribution of the center of the particle 
r is obtained by imposing that |r — rAr_i| = Rq: 



P(r) = GG(ro,r;iV) 



of Fig. S2 z 



-2a) 



Pg{p) = 27rp 



2n,ax(p) 

R-a 
(P) 



P(v/p2 + z2) dz 



R-c 



27rp 



P(Vp' + (z + 2a)2) dz 



(Sll) 



(p)+2a 



P{y/p^ + Z^) dz 



R+a 



(S12) 



as plotted in Fig. S3 3. The analytical expression for 
Po{p) is rather complicated but mean values such as 



ArjT can, in principle, be computed with- 



out major difficulty. This probability distribution can 
be simplified in two limits. For small particle radii, 
R <C \/ia, equation (S12) reads 



Pg{p) « 27r3/2pape-"('''+^o) 

X {erf[^/Q(P + a)] -erf[^/a(P-a)] 

- -erf [Va(2:max(p) + 2a)] + erf [Va2,„ax(p)] } 

« 2ape-°'p\ (S13) 

where erf(x) = lo '^V ^^"^ error function. As 
expected, we recover the Gaussian behaviour when R — >■ 
and (p) = ^-K^aLji and (p^) = 4aP/3, which corre- 
sponds to equation (10c) of Ref. [5] in the limit R <C \fLa 
(with ^ = a). 

In the large particle radius limit, R 3> \fLja, we find 
P(r) ~ Y^exp[— Q;(|r| — P)2]/(P|r|) and the most prob- 
able value for |r| get s disp laced towards non-zero values 
(|r|) = R. Equation (S12) then simplifies to 



Vg{p) 



P 



= / d'^rAT^i Gchain(0, rjv^i; iV - \)6{\v - rw^i 
exp[-a(r2-Hi?2)]^ 



3/2 sinh(2aEo|r| 



aPo|r 



-i?of 
(SIO) 



exp[-a(p^+2p2_2PVp^Tp2)]. 

(S14) 

When R increases, this probability distribution gets dis- 
placed to the right and widens as seen in Fig. [SSj p. But 
due to the projection in 2D, the most probable value of 
now scales like -y/i? and 



37r 



PVaL 



(S15) 



the 



0{\) 



a function of Irl where a 



2(Af-i)(2a)^ . The constant B 
is given by the normalization condition J^P(r)d^r = 1, 
V being the accessible volume. In order to forbid poly- 
mer trajectories intersecting the z = plane, we sub- 
stract all the forbidden paths using the mirror reflection 
argument [T] (the mirror image of z = is in the notation 



which corresponds to equation (10c) of Ref. 
limit R ^ \fLa. 

The intermediate semi- flexible regime, P/i?p 
of interest in TPM experiments, is well described by the 
worm-like chain model. For instance, in the experiments 
described in the paper, DNA lengths vary between 400 
and 2080 bp, which corresponds to 2 < L/^p < 14. This 
model can be tackled analytically [3] but it should be 
kept in mind that real chains are self-avoiding and that 
the presence of the particle renders the problem even 
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100 



FIG. S3: Theoretical probability distribution p{p) for = 100 in the two following limits, (a): Rigid rod for various particle 
sizes R: a (red), L/2 (green), L (blue), and oo (purple) (from top to bottom at u = 1). Prr(u) is expressed as a function of 
the reduced particle coordinate parallel to the substrate, u = p/pmax.(b): Gaussian chain, pg(p), p in units of a and various 
particle sizes R: [pure Gaussian see equation |S13l], 10a, 50a, 100a, and 200a from left to right. 



more intricate. For instance, the effect of the excluded 
volume of the particle is to widen and shift to large p 
the Gaussian distributions shown in Fig. |S3[ 3. It then 
becomes analytically intractable for finite chains, but it 
can be tackled numerically by Monte Carlo or Brownian 
dynamics simulations. 

In the limit L ^ where £p = 147 bp is fixed, as 
for double stranded DNA, the worm-like chain model 
tends to the Gaussian limit at a more corse-grained level 
where the Gaussian unit, the Kuhn length, is no more 2a 
but 2£p. Moreover the number of independent monomers 
N becomes L/2£p and the mean-square end-to-end dis- 
tance is Ar^ = (r^) = 2Lip. Hence the formula given by 
equation ( |S13[ ) is used where the parameter a becomes 
a = 3/(4L£p). These Gaussian distributions are com- 
pared to the numerical ones for i? = in Fig. S4 a, below. 
When i? > 0, as in Fig. |S4| 3, the numerical distributions 
cannot be compared to the present theoretical ones be- 
cause of the excluded volume effect mentioned just above. 



Equilibrium numerical distributions 

We report the equilibrium distributions computed us- 
ing Monte Carlo simulations as presented in the Material 
and Methods. 

First of all, the criterion a ^ is satisfied for all DNA 
lengths if = 50. However for small DNA molecules 
(L < 798 bp), some calculations were performed choos- 
ing = 25 to save computational time, and in this case 
we have checked that the distributions coincide within 
statistical noise. In particular, the fact that monomer 
spheres are in general larger than the actual radius of the 
DNA molecule (a > 1 nm) does not affect significantly 
our physical observables. This observation is related to 
the questions of polymer self-avoidance and solvent qual- 
ity: in our simulations, the short range repulsion, which 
is essentially of electrostatic origin, is taken into account 



by the sphere hard-core repulsion, the effective range of 
which is thus larger than the bare DNA diameter. More- 
over, there exist attractive interactions, mediated by the 
solvent, between two DNA segments which are far away 
along the chain. They are of hydrophobic nature [U [5] 
since water is a bad solvent for DNA. We assume mono- 
valent counter-ions and do not consider situations where 
DNA condensation occurs [6] . These interactions are rel- 
evant for long DNAs which are flexible at large length 
scales and can be neglected in this work, because the av- 
erage number of sphere contacts is extremely small as 
L/£p remains moderate. Therefore the choice of a ap- 
pears to be an irrelevant issue as far as chain statistics 
are concerned, provided that a <^ £p. 



The numerical distributions p{p) are plotted in Fig. S4 
Fig. S4i shows p{p) when the particle radius R = 0. Also 
displayed in the figure are the theoretical Gaussian dis- 
tributions for i? = 0, as discussed at the end of the previ- 
ous section. When L is large the numerical distributions 
and the Gaussian ones are superimposed. This observa- 
tion supports our previous remarks about the rareness 
of DNA self-contacts, as self-avoidance is not taken into 
account in the Gaussian approximation. 

Fig. |S4| 3 shows the same histograms for L — 401 and 
2080 bp where R varies between and 150 nm. For 
L — 2080 bp, the distributions look like Gaussian for all 
the values of R studied, but are slightly shifted to larger 
p and widened, as compared to the phantom chain case, 
due to self-avoidance between the particle and the chain. 
The Gaussian character is discussed afterwards. When 
L decreases, numerical distributions deviate from Gaus- 
sians: they are skewed to the right with a steep decrease 
near /Omax, a signature of the rigid-rod distribution (see 
Fig. S3 1). For L = 401 bp, this decrease is steep for small 



R, but gets smoothed when R increases. 

To quantify further the Gaussian character of the nu- 
merical distributions, we focus on their two moments (p^) 
and (p^). For a Gaussian, we expect, in 2D, (p"*) = 
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FIG. S4; Numerical distributions p{p), p in units of a [N = 50). 
(a) Symbols represent the normalized histograms for R = and 
L = 401, 798, 1500 and 2080 bp and continuous lines the theoretical 
Gaussian distributions for the same values of L. (b) Normalized 
histograms for L = 401 and 2080 bp and different values of R. Here 
lines are guides for eyes. 



2(/9^)'^. Therefore we define 



A 



,4\l/2 



(S16) 



Deviations of A from 1 are indicative of the non-Gaussian 
character. Some numerical values of A are given in the 
following table (DNA lengths L in bp, particle radii R in 
nm) . As expected, the closest values to 1 are obtained for 
long DNAs and small beads, and the Gaussian character 
is lost for short DNAs and large beads. 





L = 401 


L = 798 


L = 1500 


L = 2080 


R = 


0.85 


0.89 


0.92 


0.94 


R = 20 


0.85 


0.89 


0.93 


0.93 


R = 80 


0.85 


0.88 


0.90 


0.92 


R = 200 


0.84 


0.87 


0.87 


0.89 



Finally, in order to quantify the effect of the surface, 
we have also simulated DNAs with no surface, rotating 



freely in the whole space around their attachment point 
(data not shown). For the four DNA lengths considered 
in this work, no effect on the distribution of p has been 
detected (by contrast, z is obviously perturbed by the 
wall). Note that in both the rigid and Gaussian regimes, 
this result is shown analytically in equations (S9), (S13) 



and (S14). It appears to be also valid in the intermediate 



semi-flexible regime. 



Negligible effects of the substrate-DNA joint 

Here we address the question of the modeling of the 
substrate-DNA joint. As already discussed in Ref. [Tj in 
the different context of optically trapped particles, the 
elastic properties of the polymer's anchor point cannot 
be well characterized at the experimental level. For in- 
stance, it might happen that, for special types of exper- 
imental constructions, the joint is better modeled as a 
rigid anchor (clamped polymer). 

Therefore we examined two extreme cases: in the pa- 
per, we have considered a freely rotating joint (or freely 
flexible pivot) where all orientations of the chain tan- 
gent vector close to the substrate are allowed (except, 
of course, those intersecting the substrate). We also 
simulated a clamped joint where this tangent vector is 
prescribed to be perpendicular to the substrate by pre- 
venting the second sphere from moving (ri = 2aez). A 
prioiri, only small DNAs are likely to be affected by 
this boundary condition, because its "memory" along the 
chain is lost after a few persistence lengths. To clarify 
this question in the finite L cases considered here, we 
have also performed DMC simulations with a clamped 
joint. On the 6 parameter sets studied (namely L = 401, 
798 and 2080 bp, and i? = 20 and 150 nm), no clear ten- 
dency can be identified: within error bars, (and (r|)) 
are equal in both cases (data not shown) . Thus the exact 
way the joint is modeled is not a relevant issue. 

The observed fact that for small particles, the replace- 
ment of Dq by D± and D^^ does not modify ry, can also 
be interpreted from this perspective. In this case, the 
dynamics is dominated by the chain and weakly slow- 
ing down the particle does not change this fact. The 
only effect of hydrodynamic corrections near the wall is 
to slow down a few monomer spheres close to the sub- 
strate. This can be seen as an intermediate case between 
the two extreme previous cases: these few spheres are 
fully mobile when the joint is freely rotating, and immo- 
bile when it is clamped. Since these two extreme cases 
lead to the same relaxation times, the intermediate one - 
slowed-down spheres - also does. 



Comparison between experimental and numerical 
distributions 

Figures |S5| and |S6| show experimental and numerical 
probability distributions p(p), in the semi-flexible regime. 
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for DNA of various lengths (from 401 to 2080 bp) and var- 
ious particle radii (from 20 to 150 nm). The agreement 
is systematically good. To quantify further the small dis- 
crepancies, the values of Ary are compared into deeper 
detail in the main text. 



easily compute the relaxation time of the particle: 




100 200 

P (nm) 



FIG. S5: Experimental (black) and numerical (blue) distri- 
butions p{p) for L = 401 bp and R = 150 nm. 



Effects of exposure times Tex on experimental 
distributions 

Beyond their width Ary , probability distributions p{p) 
are also affected when Tex becomes of the order of mag- 
nitude of Til , as shown in Fig. S7 In particular, because 



of the central-limit theorem, their Gaussian character is 
restored when Tox ^ ''"y i even though the real distribu- 
tions are far from Gaussian in the rigid or semi-flexible 
regime. 

But recovering real distributions from their measured 
counterparts is a difficult task that is out of scope of the 
present work: it requires to de-convolute the averaging 
effect, which depends not only on p but also on z, be- 
cause diffusion time-scales are r-dependent. The exam- 
ination of equations (7,8) reveals that lowering the ex- 
posure time Tex is sufficient to reduce averaging effects. 
When Tox ^ Ty, they become negligible. But we have 
shown in this work that the a priori determination of the 
adapted values of Tex is not an easy task because the 
quantification of blurring effects cannot rely on simple 
assumptions. 



Simple calculation of relaxation times 

Some authors [5] estimate relaxation times ry as fol- 
lows. If one considers the complex as a spherical particle 
attached to an ideal spring (the DNA molecule), ignor- 
ing hydrodynamic corrections, excluded-volume effects, 
and interaction between the beads and the wall, one can 



2D 



(S17) 



part 



where i'part = kBT/{&TTriR) and (ry) is calculated as 
follows. If K is the spring constant of the DNA molecule, 
then 



2 fpL ' 



(S18) 



assuming a 3D phantom semi-flexible chain of persistence 
length l-p and L ^ £p 9 . Owing to the equipartition 
theorem. 



Thus one gets the theoretical value 



th.l 



47r?7^p 
knT 



RL, 



(S19) 



(S20) 



which is equal to the ratio of the friction coefficient of 
the particle, GnrjR, to the spring constant, K. 

Figure S8 (Left) shows our numerical relaxation times 



in function of RL, together with the values predicted by 
equation ( |S20 ). Note that in this case, we only compare 
our numerical values with no hydrodynamic corrections 
near the wall to this theoretical prediction, since they are 
not included in the simple theory. The above simple ap- 
proach appears to be reasonable when R is large. More 
quantitatively, if i? > L/6, then the numerical relaxation 
time is less than three times the above theoretical esti- 
mate. In addition, if R is large, the scaling Ty*^'^ ~ L is 
correct. In the best cases however, equation (S20) still 
underestimates Ty™ by a factor 2. 

By contrast, when R is small - the case of interest 
in the present work - the above approximations fail, by 
up to a factor 10. Even the scaling in equation ( S20 ) is 
incorrect since Ty™ oc in this case, because the present 
approach neglects the own dynamics of the polymer that 
dominates at small R {R < L/6). 

To improve the previous approximation, one can thus 
replace the diffusion coefficient of the particle by that 
of the particle-DNA complex, using equation (11) of the 
main text, which leads to 



th.2 



47r7;^P 



{R + L/2)L= 1 



L 
2R 



th.l 



(S21) 



in the Rouse approximation. Figure S8 (Right) shows 
that this theoretical value is better than the previous 
one since it only fails to predict the numerical values by 
a factor (at most) 2.5. The scaling is correct at small R 
since in this case Ty ^ L^. But it is now incorrect at 
larger R. 

Thus one cannot ignore the interplay between the 
(finite-size) polymer, the particle and the wall when in- 
ferring accurately relaxation times. 



17 






200 300 400 500 

P (nm) 



200 400 

P (nm) 




P (nm) 



200 400 

P (nm) 



FIG. S6: Experimental (black) and numerical (blue) distributions p(p). From top to bottom: L = 798, 1500, 2080 bp. From 
left to right: R = 20, 100, 150 nm. 



Acquisition periods Tkc 



It follows that the relative error on the measure of Tm is 



A7 



In the last section of the paper, we explain why 
the acquisition period of the detector, Tac cannot be 
much larger than the measured relaxation time, Tm. 
Indeed, if we focus for example on the correlation 
function C{t) used to measure r„i, it reads C{t) — 
C(0) exp(— t/Tm), but the only experimentally accessible 
values are C(Tac), C(2Tac), . . . Here we are looking for 
the largest values of Tac that can be used in practice. 
In this case, C(2Tac) is so close to that it does bear 
any useful information due to statistical noise. Thus the 
only relevant information, if any, comes from C(Tac) = 
C(0)exp(-Tac/T,„). If we define Q = C{T,,)/C{Q), then 
the fitting procedure essentially amounts to writing 



InQ 



(S22) 



Tm AQ 1 
Tac Q \/n' 



(S23) 



where AQ is the experimental dispersion of measures of 
Q on different trajectories and n is the number of tra- 
jectories for given experimental conditions. If Tac grows, 
then Q goes to and this relative error diverges. Thus 
Tac must not exceed a limiting value. Our less favor- 
able experimental case if for L = 798 bp, i? = 20 nm, 
and Tac = 8.9 ms, in which case Tm = 4.4 ms and 
Tm/Tac — 0.5. We have measured that AQ/Q — 0.7 
in this case. With n — 34, we get a satisfactory relative 
error of 6%. If using a higher Tac, the relative error on 
Q and thus on would rapidly increase (exponentially 
with Tac), thus requiring a much larger number n of re- 
alizations. We conclude that the condition t^/T^c ?J 0.5, 
in other words Tac ^ 2rnj, is reasonable in the present 
context. 
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FIG. S7: Experimental probability distributions p(p) for two ex- 
posure times, Tex = 5 (open symbols) and 40 ms (light red and red 
symbols). For comparison, solid symbols correspond to numerical 
distributions. Blue: R = 150 nm and L = 401 bp (t|| ~ 20 ms, 
see Table 1 in the main text); Red; R = 150 nm and L = 2080 bp 
(t|| ~ 57 ms). 



Measuring the amplitude of movement: a 
comprehensive example 

At the end of the Resuhs and Discussion section of the 
paper, we explain how the relative error on the measure 



of the amplitude of movement Ar| 
the averaging-interval length, Tav ■ 



A(Tav) 



Err(rjj) 



(rn) decreases with 



(S24) 



Figure S9 displays A(rav) measured for two experimen- 
tal data sets and two acquisition periods, Tac = 13.46 
and 40.38 ms. The second acquisition period is simply 
mimicked by discarding 2 points out of 3 in the data sets. 
The absolute error Err(r|) is calculated as the root-mean- 
square deviation of the measure of (r|) on successive in- 
tervals of duration T^v One observes that the slope —1/2 
at long times is consistent with the expected behavior. At 
short times, this behavior is potentially affected by the 
correlations between successive measurements of r|. 

More quantitatively, equation (12) of the paper pre- 
dicts an upper bound for the relative error. In the 
present case, r|| ~ 8.3 and 36.3 ms for R — 20 and 
100 nm, and T|| 2 ~ ^11/3. Thus in the case R — 100 nm 
and Tac — 13.46 ms, T|| 2 

A(Tav = Is) < ^Tjl/Tav = 

mentally A(Tav = Is) = 0.15, in agreement with this 
upper bound. In the three remaining cases, we are in the 
regime Ty 2 ^ Tac and the upper bound becomes A(Tav = 
Is) < x/Tac/Tav because successive images are uncorre- 
cted. At Tac = 40 ms and Tav = 1 s, A(Tav = Is) < 0.20, 
while we measure A(rav — Is) — 0.18 for both radii. Fi- 
nally, for i? = 20 nm and Tac — 13.46 ms, the upper 



> Tac and one expects 
0.19. We obtain experi- 




R=150 nm, no HI 
R=80 nm, no HI 
R=40 nm, no HI 
R=20 nm, no HI 
Eq. (S17) 
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FIG. S8: Left: Numerical values of the relaxation time ry™ 
in function of RL for R > (same data as in Fig. 3). The 



dotted line shows the values r| I ' predicted by equation ( S20 1 . 
Right: same data in function of {R + L/2)L. The d otted line 
shows the values Ty*^'^ predicted by equation (S21l. Log-log 
coordinates. 



bound is 0.12 and we measure 0.11. In all four cases, 
upper bound matches experimental data. 



the 



Detecting conformational changes: 
binding/unbinding 



protein 



We consider the same particle-DNA complex as in the 
example described in the Results and Discussion section 
(L = 798 bp, i? = 20 nm, = 25 and ry ~ 3 ms). 
However, instead of studying looping events, we focus on 
DNA kinking induced by protein binding on a specific 
locus. Indeed, some proteins are known to strongly bend 
the molecule when adsorbed on it, in such a way that 
binding can be detected by TPM [TUl E] . We suppose 
that the locus is situated on bead B12, in the middle 
of the chain, and that the associated bending energy is 
now ^(1 — cos(0i2 — ^kink)) instead of ^(1 — cos0i2). 
Skink is chosen equal 7r/2 in the present example. We have 
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FIG. S9: Relative error A(rav) on the measure of the ampHtude of movement (r|), extracted from experimental data for 
L = 2080 bp and ii = 20 nm (left) and 100 nm (right). Upper curves: Tac = 40.38 ms; Lower curves: Tac = 13.46 ms. Log-log 
coordinates. 



measured that in the kinked state, the particle amplitude 
of movement Ar| is reduced by about 13%: Ar| = 
0.87 Ar| ^^jjj^ijj^^j where Ary^unkinkcd = 124 nm. 

The procedure is the same as in the DNA looping ex- 
ample described in the main text. But in the present 
case, the two amplitudes are very close, which makes de- 
tection of binding/unbinding events more delicate. We 
set the threshold inbetween Ary and Ar| u,-!]^;^]^^,^ , 

at 0.93 Ary unj^iuijed- Thus relative amplitude of fluctua- 
tions must not be larger than A = 0.07 in both states to 
detect them unambiguously. Using again equation (12), 
this sets T*^ = '''II /'^^ ~ ^-^^ ^- Even though the particle- 
DNA complex is the same as in the previous example, 
one must use a much larger averaging time because the 
two state amplitudes are much closer. Again, as exampli- 
fied in Fig. SIO (left), if Tav > T*^, then false detections 



are scarce. If Tav ^ T*^, they occur with a vanishing 
probability. Note that T*^ was determined on the basis 
of a 68% confidence interval (see main text). An averag- 
ing time value four times as large as T*^ ensures a 95% 
confidence interval [equation (12)]. 

We also emphasize that in the present case, due to the 
closeness mentionned above, the two states are not distin- 
guishable on histograms p{p)- As illustrated in Fig. SIO 
(right), p{p) is not bimodal and is very similar to the 
same distribution without kinking. Whereas identifying 
bimodal distributions is useful in the context of TPM to 
determine conformational changes [10] , it is useless in the 
present context. By contrast, plotting Ar| ^^^^^^^^ vs time 
is efficient, provided that one knows in advance what av- 
eraging time Tav must be used to distinguish between 
both conformations with a good confidence. 
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FIG. SIO; Left: Simulated amplitude of movement Ary (f) for 
a L = 798 bp DNA and a 7? = 20 nm particle. The differ- 
ent plots represent Ary averaged on an interval of duration 
Tav = 0.6 s (yellow), 1.2 s (red) and 2.4 s (black). The vertical 
lines indicate the protein binding and unbinding events in the 
simulation. The horizontal solid lines show the expected val- 
ues of Ary in the kinked (bottom) and unkinked (top) states. 
The dashed line is the threshold separating the two states for 
detection purposes. Right: Distribution p{p) corresponding 
to the same simulation (solid line), together with p(p) for the 
unkinked case (dashed line, same as Fig. S4). 



